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ABSTRACT 

We present results from ASCA and ROSAT HRI observations of the powerful radio 
galaxy 4C+55.16 at redshift 0.24. Extended soft X-ray emission is imaged by the 
ROSAT HRI. The X-ray brightness profile is sharply peaked on the radio galaxy, 
characteristic of a strong cooling flow. The X-ray spectrum obtained from ASCA is 
consistent with multi-phase intracluster gas. There is evidence for an absorbed cool 
component as well as ambient cluster medium in the ASCA spectrum. A spectral fit. 
taking a cooling-flow component into account, gives a temperature of kT = 5.4lj!);g 
keV, metal abundance (0.5 ± 0.1)^0, excess absorption on the cool component, 
ATVh — 4.9^^3 X lO^^cm^^, and an absorption-corrected bolometric luminosity of 
2.2 X 10'*^erg s^^ (_ffo = 50 km s^^ Mpc^^, go = 0.5). The mass deposition rate is 
estimated to be ~ llOOM© yr^^ from the spectral analysis, in good agreement with 
that derived from imaging analysis of the ROSAT data when corrected for absorp- 
tion. We tentatively identify a blue feature, seen ~ 15 arcsec from the centre of the 
radio galaxy in a published optical image, as a gravitationally-lensed arc. The inferred 
lensing mass is consistent with the gravitational mass derived from the X-ray data. 
The best-estimate of the redshift of the lensed object is ^ 1.5(> 0.7). All the observed 
properties suggest that the environment of 4C+55.16 is similar to known massive cool- 
ing flow clusters. This is the first massive cooling flow to be found around a powerful, 
radio source with a compact, GHz-peaked spectrum core. 

Key words: galaxies: clusters: individual: 4C-I-55.16 - X-rays: galaxies - cooling 
flows - gravitational lensing 



1 INTRODUCTION 

4C-I-55.16 (= 0831-1-557) is a radio galaxy at a redshift of 
0.240 (Lawrence et al 1996) in the Pearson- Readhead (1981, 
1988) 5 GHz flux-density-limited complete sample. The ra- 
dio source is compact and powerful (1.1 x 10^^WHz~'^sr~^ 
at 5 GHz). In the high resolution radio maps obtained 
from MERLIN and VLBI (Whyborn et al 1985; Pearson 
& Readhead 1988), highly irregular, small-scale structures 
have been resolved. The radio spectrum of the compact core 
shows a turnover around 1 GHz (Component C in Why- 
born et al 1985), resembling the class of Giga-hertz Peaked 
Spectrum (GPS) sources (e.g., O'Dea 1998) whilst the mini 
double-lobe with a projected size of 11" (54 kpc at a red- 
shift of 0.24) has a steep, low-frequency spectrum. A sum 
of the two components results in a flat radio spectrum over 
the GHz range. 

A blue continuum contributes ~ 50 per cent of the light 
at 3850A (Heckman et al 1983) and its optical emission-line 



spectrum shows a medium ionization state (e.g., Heckman et 
al 1983 based on the [OII]A5007/[OII]A3727 ratio; Whyborn 
et al 1985; [OIII]A5007/H/3 ~ 4, Lawrence et al 1996). 

X-ray emission was detected during the ROSAT All Sky 
Survey (RX J08349-I-5534, Brinkmann et al 1995; Laurent- 
Muehlei et al 1996; Bade et al 1998). The ROSAT papers 
assumed that the X-ray emission originates in an active nu- 
cleus residing in 4C-I-55.16. However, as our imaging and 
spectral study using the ROSAT HRI and ASCA data show 
below, the observed X-rays appear to be dominated by clus- 
ter emission surrounding the radio galaxy. 

The optical CCD images of 4C-f55.16 were obtained 
by Hutchings, Johnson & Pyke (1988) using the Canada- 
France-Hawaii Telescope (CFHT) with B and R filters. 
There is a blue feature about 15 arcsec southwest from the 
centre of the galaxy, which we tentatively identify as a gravi- 
tationally lensed arc. We find a good agreement between the 
cluster masses estimated from X-ray and lensing techniques. 
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Figure 1. The X-ray (0.1-2.4 kcV) contours of 40+55.16 ob- 
tained from the ROSAT HRI image overlaied onto the optical 
DSS image talcen by Palomer 48-inch Oschin Schmidt Telescope. 
The DSS image is not deep enough to show the arc-like feature de- 
tected in the CFHT image (Hutchings et al 1988). Contour levels 
are 6, 10, 16, 25, 40, 63, 95 per cent of the peak brightness which 
is 63 ct pixel^^. The estimated average background is about 1.3 
ct pixel" ^. 



of the HRI. The raw HRI image has been smoothed with 
the adaptive kernel method, ASMOOTH (Ebeling, White 
& Rangarajan 1998), using a gaussian kernel and a charac- 
teristic smoothing threshold, above the local background, of 
4(7. Fig. 1 sows the X-ray intensity contours of the smoothed 
data overlaid on the optical Digitized Sky Survey (DSS) im- 
age. 

The X-ray emission is sharply peaked at the centre of 
the extended image, as observed in strong cooling flow clus- 
ters. The X-ray peak is displaced from the galaxy nucleus 
by ^ 7 arcsec, well within the mean displacement (25 arc- 
sec) between X-ray and optical peaks seen in cooling flow 
clusters observed with ROSAT HRI (Peres et al 1997). Al- 
though the displacement is consistent with the positional 
uncertainty of ROSAT pointing (^ 10 arcsec), if it is real 
the galaxy resides in the X-ray cavity, as seen in the Perseus 
cluster (Bohringer et al 1993). There are three point-like 
sources detected in the HRI field. One of them shows an off- 
set from a possible optical counterpart, similar in amplitude 
and direction to that for 4C-f55.16. Therefore the displace- 
ment between X-ray and optical peaks may be merely due 
to a pointing error. 

The elongation of the X-ray image is prominent at the 
low suface brightness levels. However, this may be partly 
due to a contamination from point sources in the field, par- 
ticularly in the SW region. The present observation is too 
short to investigate such a faint X-ray morphology. 



2 OBSERVATIONS AND DATA REDUCTION 

404-55.16 was observed with ASCA (Tanaka, Inoue & Holt 
1994) and ROSAT HRI (Pfeffermann et al 1987). ASCA pro- 
vides moderate resolution imaging and spectra in the 0.5-10 
keV band while the ROSAT HRI provides high resolution 
(FWHM ~ 5 arcsec) 0.1-2.4 keV imaging. A summary of 
the observations is given in Table 1. 

The four detectors onboard ASCA, the Solid state 
Imaging Spectrometers (SIS; SO and SI), the Cas Imaging 
Spectrometers (CIS; G2 and G3) were operating normally. 
The radio galaxy was pointed at the nominal position on the 
best-calibrated CCD chip of each SIS detector. The lowest 
energy event threshould in the SIS was set at 0.47 keV. Stan- 
dard calibration (used for the Rev2 processing) and data re- 
duction techniques were employed, using FTOOLS (version 
4.0) provided by the ASCA Cuest Observer Facility at Cod- 
dard Space Flight Center. The SIS data were corrected for 
the effects of the Dark Frame Error (DFE) and 'Echo' and 
hot/flickering pixels were removed. 

The ROSAT HRI observation was short (16.6 ks). The 
soft X-ray image revealed extended X-ray emission which 
could not be resolved with ASCA (Section 3). The data 
were also analysed by the deprojection technique (Section 
5). There is no evidence for X-ray variation in both ASCA 
and ROSAT data. 



3 THE ROSAT HRI IMAGE 

The X-ray source detected at the position of 4C+55.16 is 
clearly extended beyond the Point Spread Function (PSF) 



4 THE ASCA X-RAY SPECTRUM 

The spectral analysis was performed using XSPEC (version 
10.0). The MEKAL model (the original MEKA code, de- 
scribed by Kaastra 1992, with modified Fe-L line emissivity 
by Liedahl et al 1995) for an optically thin, coUisional ion- 
ization equilibrium plasma was used with solar abundances 
taken from Anders & Grevesse (1989). The photoelectric 
absorption model was taken from Morrison & McCammon 
(1983). The Galactic absorption at the position of 4C-I-55.16 
is estimated to be A^h ~ 4.2 x 10^" cm^'^ from the HI mea- 
surements of Dickey & Lockman (1990). Absorption column 
densities obtained from the spectral fits are then excesses 
above the Galactic value. Quoted errors to the best-fit spec- 
tral parameters are 90 per cent confidence regions for one 
parameter of interest. 

The ASCA spectrum shows a clear line feature at 5.4 
keV, which is in an excellent agreement with the rest Fe Ka 
line emission at ~ 6.7 keV expected from a thin thermal 
plasma with a temperature of several keV at the redshift 
{z = 0.240) of 4C+55.16. The agreement between the red- 
shifts of the Fe Via line and the radio galaxy strongly sup- 
ports the presence of a cluster around the radio galaxy rather 
than in the background which was suspected by Hutchings et 
al (1988). The thermal emission model (MEKAL) provides 
a slightly better fit to the data (0.6-9 keV from the SIS; 
0.9-10 keV from the GIS) than the model of a power-law 
plus a gaussian line for the Fe Ka line feature (see Table 2). 
The strong, narrow Fe Ka line at 6.7 keV is unlikely for an 
active galaxy but naturally explained by thermal emission 
from cluster gas whose X-ray emission has been spatially 
resolved by the ROSAT HRI. 

There is evidence for multi-phase gas in the ASCA 
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Table 1. The ASCA and ROSAT observations of 4C+55.16. No eorrections for vignetting has been made for the count rates given here. 



Satellite 


Detector 


Operation mode 


Band 


ASCA 


SIS (SO/SI) 


ICCD/Faint 


0.5-10 keV 




CIS (G2/G3) 


PH 


0.7-10 keV 


ROSAT 


HRI 




0.1-2.4 keV 



Date 



Exposure 
ks 



Count rate 
cts^i 



1996 Nov 1-2 34.5/34.0 0.12/0.10 
38.1/37.8 0.065/0.078 

1997 Oct 20 16.5 0.062 



Table 2. Spectral fits to the ASCA SIS (0.6-9 keV) and CIS (0.9-10 keV) data. ATVh is the excess column density above the Galactic 
value (4.2 X lO^'^cm"^), measured in the galaxy rest frame. Errors are quoted at the 90 per cent confidence level for one parameter of 
interest. In the model (a), the gaussian line energy is corrected for the galaxy redshift. The equivalent width of the line is 520_j3g eV. 
In the thermal emission model (MEKAL), the Solar abundance table from Anders & Grevesse (1989) is used. In the model (c), the cool 
component is the cooling flow model with an initial temperature kTi, metal abundane Z and a mass deposition rate of M. The column 
density of cold absorption for the cool component (ATVjj cool) is measured a,t z = 0.24. 
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(a) Power-law plus a Gaussian line 
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(b) Single thermal emission model (MEKAL) 



Z 

Zq 



0.48 



ANn 
lO^^cm^^ 
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+0.11 
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(c) Multi-phase thermal emission model 



0.49 



+0.13 



M 
Meyr-i 
11QQ+240 



AAfH.cool 

102lcm-2 



4.9 



+3.4 
1.3 



xVdof 
529.6/511 

xVdof 
521.8/513 

xVdof 
524.0/512 



spectrum. Fits to the ASCA data above and below 3 keV 
with a single thermal emission model (MEKAL) give sig- 
nificantly diflerent values of temperature (Table 3). The 
data below 3 keV also require excess absorption (AA''h ~ 
0.98+032 ^ lO'^^cm"'^ measured in the Earth frame) above 
the Galactic value. The extrapolation of the best-fit model 
for the 3-10 keV data leaves excess emission down to 1 keV 
followed by a decline towards 0.6 keV is probably due to ab- 
sorption (Fig. 2) . This indicates a multi-phase gas consisting 
of at least an absorbed, cool component with a less absorbed, 
ambient medium. This is characteristic of a cooling flow. 



A fit to the whole band data (0.6-8 keV for the SIS; 
0.9-10 keV for the GIS) with a muhi-phase model (Model- 
c in Table 2), consisting of a single temperature MEKAL 
with Galactic absorption and the cooling flow model (John- 
stone et al 1992) modifled by extra absorption at the source, 
gives temperature, metal abundance, mass deposition rate 
and excess absorption column for the cooled component of 



kT = 5.4^ 



keV, Z 



:O.5Oto:l3^0,M 



1100 



+240 



Meyr- 



and A'^H ~ 4.9+'J'3 x lO^^cm"^, respectively. The temperature 
is listed for the ambient medium (i.e., before cooling), the 
metallicity is assumed to be identical between the two com- 
ponents, and the absorption column density is corrected for 
the galaxy redshift. The quality of the fit is acceptable and 
comparable to the single phase model (see Table 2). The 




Energy (keV) 

Figure 2. The ratio plot of the data and the extrapolation of the 
MEKAL model best-fitting the 3-10 keV data. Excess soft X-ray 
emission can be seen down to 1 keV and then it declines towards 
lower energies. 



covering fraction of the cold absorption must be larger than 
0.9 (90 per cent lower limit). 

The observed fluxes in the 0.5-2 keV and 2-10 kcV 
bands obtained from the GIS are 1.72 x lO^^^erg cm^^ s^^ 
and 2.56 x lO^^'^erg cm~^ s~^, respectively. In the best-fit 
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Table 3. Thermal emission model fits to the data above and below 3 keV. 



Data 

E>3 keV 
E<3 keV 



kT 
keV 



5.40 
3.73 



+ 1.01 
0.81 

+0.71 
0.46 



0.48 

0.59 



=-0 

+0.12 
0.11 

+0.37 
0.25 



lO^^cm" 



0.98 



+0.34 



x'/dof 

151.4/156 
406.7/394 



ASCA SIS 
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Figure 3. The Fe K band spectrum obtained from the ASCA SIS. 
The solid-line histrogram represents the best-fit thermal emission 
model (See Table 3). Fe K/3 emission around 6.4 keV (7.9 keV 
in the galaxy rest-frame) is detected with marginal significance 
above the model. 



multi-phase model, about half of the total flux comes from 
the cool component. 

Excess line-like emission at a rest energy of 7.9 keV 
(6.4 keV in the observed frame) is marginally detected in 
the SIS data (Fig. 3). It can be identified with Fe K/3. This 
feature is barely seen in the G3 data but not in the G2 data 
which show an unusual deficit between 6-7 keV, probably 
due to some anomaly in the detector. Since there is an in- 
strumental feature at 6.4 keV in the SIS spectrum, the SIS 
results should be treated with caution. If the feature in the 
SIS data is real, the Fe K/3 emission is underestimated by a 
factor of 2.3 ± 1.5 by the best-fit thermal emission model. 
An anomalous Fe K/3/Fe Ka ratio has been observed in the 
central region of a few cooling fiow clusters, and has been 
interpreted as the effect of resonant scattering of line pho- 
tons in the cluster core (e.g., Akimoto et al 1996; Molendi 
et al 1998). This interpretation is however not applicable to 
4C-I-55.16 because the whole cluster emission is observed. 
Further deeper observations are required to confirm the ex- 
istence of this emission feature. 



5 DEPROJECTION ANALYSIS OF THE HRI 
DATA AND THE X-RAY MASS MODEL 

We have carried out a deprojection analysis of the ROSAT 
HRI data. An azimuthally-averagod X-ray surface brightness 



profile was constructed for the cluster. This was background- 
subtracted, corrected for telescope vignetting and re-binned 
into 12 arcsec bins to provide sufficient counts in each annu- 
lus for a reliable statistical analysis to be carried out. With 
the X-ray surface brightness profile as the primary input, 
and under assumptions of spherical symmetry and hydro- 
static equilibrium, the deprojection technique yields the ba- 
sic properties of the intracluster gas (temperature, density, 
pressure, cooling rate) as a function of radius. 

The deprojection method requires the total mass pro- 
file for the cluster to be specified. We have iteratively de- 
termined the mass profile (which has been parameterized as 
an isothermal sphere; Equation 4-125 of Binney & Tremaine 
1987) that results in a deprojected temperature profile which 
is isothermal within the region probed by the ROSAT data 
and which is consistent with the best-fit temperatures deter- 
mined from the ASCA spectra, using the cooling-fiow model 
(Section 4, The validity of the assumption of isothermal 
mass-weighted temperature profiles in the cluster cores is 
discussed by Allen 1998b). The best-fitting mass model has 
a core radius of 60 ± 20 kpc and a velocity dispersion of 



820: 



km s 



The primary results from the deprojection analysis 
are as follows: for an assumed Galactic column density of 
4.2 X 10^"cm~^, we determine the mean cooling time within 
the central 12 arcsec bin of tcooi = 2.0lg'2 ^ ^0^ Y^^ ^ cool- 
ing radius (beyond which the cooling time exceeds a Hubble 
time) of rcooi = 180140" kpc, and an integrated mass de- 
position rate within the cooling radius of M = 46O1^4qM0 
yr~^. If we correct for intrinsic absorption in the cluster, 
as determined from the ASCA data, these values are ad- 



justed to icool 
M 



1.5t?,i X 10^ 



yr, rcooi = 270li(, kpc, and 



970+11°Mq yr-\ 



6 LENSING ANALYSIS AND COMPARISON 
WITH THE X-RAY RESULTS 

The mass model determined from the X-ray data may be 
compared to the mass implied by the observed lensing con- 
figuration in the cluster (Section 1). Since only a single, 
putative gravitational arc is seen, and to be consistent 
with the X-ray analysis, we have only carried out a sim- 
ple, spherically-symmetric analysis of the lensing data. For 
a spherical mass distribution, the projected mass within the 
tangential critical radius, which we assume to be equal to 
the arc radius, rare = 15 arcsec (71.8 kpc), is given by 



.if^arc l^^arc J 



Da 



4:(_T \ -Lucius -L-'arc — clue 



(1) 



where Dcius, i^arc and Darc-cius are respectively the angu- 
lar diameter distances from the observer to the cluster, the 
observer to the lensed object, and the cluster to the lensed 
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Figure 4. The projected Icnsing mass within the critical radius 
as a function of the redshift of the lensed object (solid curve). 
The gravitational mass estimated from the deprojection anlaysis 
using the ASCA results is indicated by a dashed line with the 90 
per cent confidence limits shown by dotted lines. 



object. Fig. 4 shows the projected mass within the critical 
radius as a function of the redshift of the arc (solid curve) . 
The horizontal dashed and dotted lines mark the best fit 
(projected) mass measurement and 90 per cent confidence 
limits determined from the X-ray data, within the same ra- 
dius (3.8tog X 10^^ Mq). We see that the X-ray and lens- 
ing mass measurements are consistent for any arc redshift 
Zarc > 0.7. The best match between the X-ray and lensing 
mass measurements is obtained for an arc redshift of 1.5. 



7 DISCUSSION 

The ROSAT HRI image of the powerful radio galaxy 
4C-I-55.16 shows extended X-ray emission peaking at the ra- 
dio galaxy indicating cluster emission with a strong cooling 
flow. 

A spectral study of the ASCA data suggests the X-ray 
emitting gas to be multi-phase. An absorbed, cool compo- 
nent is found in the spectrum. The muti-phase spectral anal- 
ysis indicates that the temperature of the ambient cluster 
medium is kT ~ 5.4 keV. A single-phase model fitted to 
the data gives a temperature lower by ~ 1 keV, typical of a 
cooling flow cluster. 

The mass deposition rate of the cooling flow, 
1100^410^0 y^~^ y derived from the spectral analysis is con- 
sistent with that (97OI45OM0 yr~^) estimated from the im- 
age analysis when corrected for excess absorption. Agree- 
ments between mass deposition rates derived from the two 
methods have been found for other distant cooling-flow clus- 
ters (Allen 1998b). 

The optical spectrum of 4C+55.16 (Lawrence et al 
1996) is indeed very similar to the other cooling flow galax- 



ies (Crawford et al 1999). The Hq luminosity is about 
8 X 10^^ erg s^^ (Lawrence et al 1996). The relatively large 
Balmer decriment (Hq/H/3 ~ 5.5) suggests a significant red- 
dening, which is often observed in cooling flows. The inferred 
absorption column density is slightly smaller than observed 
in the X-ray spectrum. 

The absorption-corrected 2-10 keV and bolometric lu- 
miosities, computed from the multi-phase model, are 8.0 x 
10*''erg s"\ and 2.2 x lO^'^erg s~\ respectively {Hq = 50 
km s^^ Mpc^^ and go ~ 0.5). About 60 per cent of the 
bolometric luminosity is due to the cooling flow. The bolo- 
metric luminosity exceeds that predicted for the single-phase 
temperature of 4 keV from the correlation between emission- 
weighted cluster temperature and luminosity (Mushotzky 
1984; Edge & Stewart 1991; David et al 1993; Fabian et 
al 1994; Mushotzky & Sharf 1997; White et al 1997). A 
similar descrepancy is found for the other strong cooling 
flow clusters (e.g., Fabian et al 1994; Allen & Fabian 1998a; 
Markevich 1998). Taking the temperature derived from the 
multi-phase spectral analysis, 4C-f55.16 fits well the fcTx- 
LboI correlation obtained from a similar analysis of other lu- 
minous (I/Boi > 10''^ erg s~^) clusters for which the effect of 
cooling flows is included (Allen & Fabian 1998a), and is con- 
sistent with Z/Boi oc Tx expected from simple gravitational 
collapses for formation of clusters (Kaiser 1986; Navarro, 
Frenk& White 1995). 

As shown in Section 6 (and Fig. 4) , the mass estimated 
using the tentatively-identified lensing arc and the gravita- 
tional mass derived from the X-ray deprojection analaysis 
are in good agreement. Moreover, the core radius of 60 ± 20 
kpc measured for 4C+55.16 is similar to the best-fit mean 
value of ~ 50 kpc measured in the six lensing cooling-fiow 
clusters studied by Allen (1998b). 

The inferred metallicity of this cluster gas is ~ 0.5^©, 
which is not unusual, but is certainly one of the higher val- 
ues measured among the ASCA cluster sample compiled 
by Allen & Fabian (1998b). They showed that cooling-fiow 
clusters show higher metallicity than non cooling-fiow clus- 
ters, and suggest that the sharply peaked X-ray brightness 
profiles may give apparently high values of the emissivity- 
weighted metallicity in cooling flow clusters when there is a 
metallicity gradient in the cluster core. Thus the high metal- 
licity in 4C-I-55.16 may be due to a steep metallicity gradient 
towards the cluster centre. The high metallicity measured in 
the spectrum also rules out significant contribution from the 
active nucleus to the observed hard X-ray emission around 
the Fe K band, otherwise the line would be less prominent. 

We conclude that all of the available evidence points to 
the environment of 4C-I-55.16 being like other distant mas- 
sive cooling fiow clusters. 4C+55.16 is yet another power- 
ful radio galaxy surrounded by a strong cooling fiow. Un- 
like Cygnus A and 3C295, 4C-I-55.16 is a compact radio 
source. Detailed high spatial resolution observations with 
AXAF will be required to determine the interaction between 
the radio source and the dense surrounding intracluster 
medium. Probably the most comparable cluster to 4C-I-55.16 
is PKS 0745-19, which also contains a strong cooling fiow 
and exhibits gravitationally-lensed arcs (Allen, Fabian & 
Kneib 1996). However, the amorphous radio source shows 
a steep spectrum and is almost 2 orders of magnitude less 
powerful thtan 4C+55.16 at 5 GHz. Massive cooling fiows 
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have not been found so far around 'pure' compact, GPS 
cores (e.g., O'Dea et al 1996). 

The gas pressure at the centre of a coohng flow is P = 
nT Ri 10^ cm"^ K (Heckman et al 1989). This leads to 
a free-free absorption optical depth rg ~ O.GP^iytJ^T^ h 
at a frequency 10 ug Hz, through a cloud of length 10 Zi 
pc at pressure lO^P? cm""^ K and temperature 10*r4 K. 
It is therefore plausible that the turnover seen in the radio 
spectrum of the compact core at ~ 1.7 GHz (Whyborn et al 
1985) is due to free-free absorption in the Ha emitting gas 
close to the nucleus. The covering fraction in such clouds 
must however decrease rapidly away from the nucleus in 
order that the radio knot observed ~ 100 pc north-west of 
the nucleus (Whyborn et al 1985; Pearson & Readhead 1988) 
is not absorbed. The inverted spectrum of the counter-jet of 
NGC1275, which also resides in a cooling flow, is suspected 
to be due to free-free absorption on the pc scale (Vermeulen, 
Readhead, & Backer 1994). Free-free absorption is favoured 
for the spectral turnover in GPS sources by Begelman (1997) 
and Bicknell, Dopita & O'Dea (1997). Part of the observed 
optical emission-line luminosity from the innermost part of 
4C-I-55.16 can be expected from shocked gas surrounding the 
jets of such a powerful and a relatively young radio source 
(e.g., Bicknell et al 1997), as well as from the cooling flow. 
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